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I. INTRODUCTION 



In our previous paper [|l], we have formulated in formal and abstract ways M-fold su- 
persymmetry in quantum mechanics p|, H, ||, |[ and investigated general properties of 
the models which possess this symmetry. A/'-fold supersymmetry is characterized by the 
supercharges which are A/'-th order polynomials of momentum and similar generalizations 
of supercharges were also investigated in different contexts ||, ^, 0, |TT], 0, |13|, |T^, |T5|, |T^, 
p!7| , |l^, 19, ^ 21, 22, 2^. We have shown that A/'-fold supersymmetric models have a lot of 



significant properties similar to the ordinary supersymmetric ones [^, ^ ^ |2^ such as de- 
generate spectral structure between bosonic states and fermionic ones, non-renormalization 
theorems for the generalized Witten index and for a part of the spectra, and so on. Fur- 
thermore, we have introduced the notion of quasi- solvability to identify an important aspect 
of A/'-fold supersymmetry and have proved the equivalence between A/'-fold supersymmetry 
and quasi-solvability. Recently, we have further shown ^ that Type A subclass of A/'-fold 
supersymmetry which was first introduced in Ref. 0| is equivalent to the quasi-solvable 



models constructed by s/(2) generators |]28[. Then, it has turned out that the equivalence 
between them for special cases which had been reported previously in Refs. [@, 0, is 
generically holds. 

Quasi-solvability means the existence of a finite dimensional invariant subspace under 
the action of the Hamiltonian. As a consequence, a part of the spectra can be solved by a 
finite dimensional algebraic equation. In the case where the subspace is physical, that is, L^, 
these spectra give the true eigenvalues of the Hamiltonian. In this case, the system is often 
called quasi-exactly solvable p8|, On the other hand, if the subspace is not physical, 

solvable spectra only give perturbative eigenvalues at most and thus we have dubbed this 
case quasi-perturbatively solvable This distinction is quite important, especially in view 
of dynamical A/'-fold supersymmetry breaking; A^-fold supersymmetry is broken dynamically 
if a system is quasi-perturbatively solvable while it is not broken if a system is quasi-exactly 
solvable. 

In this article, we analyze particular models which can be A/'-fold supersymmetric more 
quantitatively. Previously in Ref. 0], an asymmetric quartic double- well potential was 
investigated in detail by the valley method [§T[ |3|, |S|, It was 

shown that the system possesses A^-fold supersymmetry at specific values of a parameter 
in the Hamiltonian where the leading Borel singularity of the perturbative corrections for 
the first A/'-th energies disappear. This result is a consequence of a general property of A/'- 
fold supersymmetry, namely, the non-renormalization theorem. It was also shown that the 
non-perturbative corrections for the first A/'-th energies do not vanish even when the system 
becomes A/'-fold supersymmetric. This result consistently reflects the fact that, in the case 
of an asymmetric quartic double-well potential, the solvable subspace is not physical, that 
is, the system is quasi-perturbatively solvable; non-vanishing non-perturbative effects break 
A/'-fold supersjTiimetry dynamically. 

These observations show that combining with the formal discussions in Ref. |l|] quanti- 
tative analyses may give deeper and complement understanding on dynamical properties of 
the A^-fold supersymmetric models. As in the case of the ordinary supersymmetric models, 
non-perturbative analyses are quite important in the A/'-fold supersymmetric case; dynami- 
cal A/'-fold supersymmetry breaking can take place via purely non-perturbative effects, e.g., 
quantum tunneling. However, non-perturbative analyses are in general quite non-trivial 
even in the simple one-dimensional quantum mechanics. The valley method is one of the 
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most successful tools for this kind of purpose. So we fully employ it in this work. 

The article organizes as follows. In the next section, we summarize the general results and 
properties of the Type A subclass of A/"- fold supersymmetry [jl], |, ^, ^. In section |T|, we 
develop particular cases of the Type A models which are especially relevant for the analyses 
in this article. Sections ^ and |V| are devoted to valley method analyses on a periodic and a 
sextic triple-well potential, respectively. We choose potentials to be investigated so that the 
systems can be Type A A/'-fold supersymmetric at specific values of a parameter involved 
in the potentials. The way of the choice enables us to clarify the characteristic features of 
Type A A/'-fold supersymmetry. The periodic potential is always quasi-exactly solvable when 
Type A A^-fold supersymmetric while the triple- well potential can be either quasi-exactly 
or quasi-perturbatively solvable. In both the cases we show that the disappearance of the 
leading Borel singularity occurs. However, the non-perturbative corrections vanish when 
and only when the systems are quasi-exactly solvable. Finally, we give summary in the last 
section. 

II. GENERAL PROPERTIES OF TYPE A AA-FOLD SUPERSYMMETRY 

First of all, we summarize the general results and properties of Type A A/'-fold super- 
symmetry. For details, e.g., derivation of the results, see Refs. |^. To define A/'-fold 
supersymmetry, we introduce the following Hamiltonian and the A/'-fold supercharges, 

H^^ = H]^{p, g)^^t + H^{p, g)^t^^ (2.1) 

QAr = P|(p,g)^, Q^j^ = PM{p.qW- (2.2) 
Here ip and ip'^ are fermionic coordinates which satisfy, 

{^,^} = {^t^^t} = 0, {V',V'^} = 1, (2.3) 
and are usually represented as the following 2x2 matrix form: 

The component of the A/'-fold supercharges Pj\j- is given by an A/'-th order polynomial of 
p = —id/dq and thus expressed as, 

Pm=P^ + wx-i{q)p^~^ + ■■■ + wi{q)p + Wo{q), (2.5) 

without any loss of generality. Then, the system is said to be A/'-fold supersymmetric if the 
Hamiltonian H^v" commutes with the A/'-fold supercharges Q^f and Q\f: 

[Q^f,U^] = [Q^^,U^] = 0. (2.6) 
The Type A A/'-fold supersymmetry is characterized by a particular class of the A/'-fold 



3 



supercharges which can be expressed as the following form:^ 

(A/'-l)/2 

= n {D + ikE{q)y D = p-iW{q). (2.7) 

fe=-(Ar-l)/2 

If we restrict to be the following Schrodinger type, 

H% = \p' + V^{ql (2.8) 

we can show [^] that necessary and sufficient conditions of the Hamiltonian ( p.l| ) with ( p.8[ ) to 
be Type A AT-fold supersymmetric, that is, to satisfy the relation ( p.6[ ) , are as the following: 



V^{q) = \w{qf + ^E{qf - 2E'{q)) ± :^W^'(g), (2.9a) 

(^^ - 2E(g)) (^^ _ i5(5)^ ^ (^^ + i5(g)^ i5(g) = (AT > 3). (2.9c) 



A. The Solvable Subspaces 

Owing to the relation Eq. (|2.6| ), the A/'-dimensional vector spaces defined by 

V^ = kerP^, V_^ = kerP| (2.10) 

are invariant under the action of Hj^ and if^, respectively. We can therefore define the 
matrices as follows: 

= J] S±^</.±, (2.11) 

m=l 

where 0^ are bases of the V^, respectively. It can be proved [|I| for general AA-fold supersym- 
metry that the mother Hamiltonian Tij^ defined by the anti-commutator of the supercharges 
can be expressed as, 

7^,, = i|nt, n,A-UdetM+,(i7^)+p+P| A , . 

UN- ^iQM^Qjsf\- detM^(/f_^)+p-tp^y'' ^^-^^^ 



where, 



M±(A) = 2(AI-S±), (2.13) 



^ Note that W{q) in this paper is the same as W{q) and not as W{q) in the old notation of our previous 
paper ^ ^, ^. Since it is a bit troublesome to keep the unnecessary tilde, we omit it. 
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and p"^ are at most (TV — l)-th order differential operators. From the definition of the 
mother Hamiltonian ( |2.12| ), the elements of the subspaces are also characterized as the 
zero-modes of the mother Hamiltonian. 

In the case of Type A, we can obtain analytic expressions for these bases: 

€{q) = hiqr-'h'iq)-^^-'y'U{q)^\ (n = 1, . . . , AT), (2.14) 

where, ^ 

U{q) = e^^''^(i\ (2.15) 
and h{q) is a solution of the following linear differential equation: 

h'\q)-E{q)h\q)=0, (2.16) 

and thus generically given by, 

h{q)=ci ! dqe^'^^^^'^'^ + C2. (2.17) 



The appearance of the two arbitrary constants Ci^2 in Eq- ( |2.17|) reflects the fact that the 
spaces spanned by the bases Eq. ( |2.14| ) are invariant under any linear transformations 



on h{q). With the aid of these bases Eq. ( |2.14| ), the components of the matrices S^^ defined 



by Eq. ( |2.11[ ) can be determined (for each fixed = 1, . . . , A/") by the following recurrence 
relations: 



P. 



for m = 1, . . . , A/" — 1 with the initial conditions. 



From Eq. ( p.ll| ), the spectra of the Hamiltonians in the subspaces are given 

by, 

detM±(i?±) = 0. (2.20) 

If (/)n(g)'s are normalizable, linear combinations of them which diagonalize the matrix S 
are the true eigenstates of Hj^. In this case, the system is often called quasi- exactly solv- 



able ||2^, |29|, 0. On the other hand, if 0„(g)'s are not normalizable, they have, at most, 
restricted meanings in the perturbation theory. In this case, the spectra determined by 
Eq. (p.20|) only give perturbatively correct ones. For this reason, we dub the case quasi- 
perturbatively solvable [Q. Then, AA-fold supersymmetry of the total system is dynam- 
ically broken when both of the systems are quasi-perturbatively solvable. Otherwise, 
that is, at least one of the systems is quasi-exactly solvable, the elements of the cor- 
responding solvable subspace give the A/'-fold supersymmetric physical states and therefore 
A^-fold supersymmetry is preserved. 



The definition of U{q) is also different from the one in our previous papers. 
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B. A Non-renormalization Theorem 



A kind of the non-renormalization theorem holds for the Type A models. We first assume 
that we can set W{0) = by the redefinition of the origin of the coordinate q and the energy. 
To define a perturbation theory, we then introduce a coupling constant g as, 

Wiq) = -wigq), E{q) = g e{gq), (2.21) 
9 

so that, in the leading order of g, the potential become harmonic with frequency |w'(0)|, 

V^{q) = \w\Qfq' + 0{g). (2.22) 



From Eq. ( ^.14| ) with Eqs. ( ^.15| ), ( ]2.17| ) and ( |2.21D , we can easily see that the 0^ behave 
as, 

0^(9) = U{Qt' (g"-i + 0(^?))e±-'(°)'^'/l (2.23) 



Here, we choose the two arbitrary constants ci,2 in Eq. ( p.l7|) as /i(0) = 0, h'iOi) = 1. Thus, 



as far as w'{0) > (< 0), all the 0~(0^) remain normalizable in any finite order of g even 
if 4>n{4'n) themselves are not normalizable. So, they stay the TV-fold supersymmetric vacua 
in any order of the perturbation theory and therefore any perturbative corrections do not 
break A/'-fold supersymmetry. 



III. SPECIAL CASES OF TYPE A 7^-FOLD SUPERSYMMETRY 

In this section, we illustrate some special cases of the Type A A/'-fold supersymmetry by 
using the general results obtained in the previous section. 

A. Exponential Type Potentials 

At first, we will consider the case where E{q) = A (non-zero constant). This is a trivial 
solution of Eq. ( p.9c|) . From Eq. (|2.9b|) we yield. 



W{q)=Cie^'^ + C2e-^'^ + Gi. (3.1) 
In this case, the Hamiltonians and the supercharge are given by 

H^=-p' + -W{qr + ^^X'±'^W'{q), P^= n {D + ^kx). (3.2) 

k=-{J\f~l)/2 

The function h{q) can be chosen as, 

h{q) = ^. (3.3) 
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Bases of the solvable subspaces are calculated as, 



X 



n-l 



exp 



l(Ar -2n + l)\q ± Csq ± y e"^ T ^ 



(3.4) 



Thus, normalizability of 0^ depends on the values of the constants Ci and A. For example, 
provided that all the constants Ci, C2 and A are non-zero real numbers, either 0"*" or (f)~ is 
normalizable when C1C2 < while both of 0^ are not when C1C2 > 0. The correspondence 
between quasi-exact solvability and A/'-fold supersymmetry in the case of the exponential 
type potentials Eq. 



|T|) was recently discussed in Ref. [p!^ . 
The non-zero matrix elements of can be calculated as follows. The direct action of 
the Hamiltonians Eq. (|2.8| ), with the Type A potentials Eq. (p.9a|) , on the bases Eq. (p.l4|) 
reads. 



1. 



l)(n-2)/i' 



{Af -2)h" T^Wh' 



1 



12 



[{Af - 2) {E' + E'^)tQ {W + EW 



From Eqs. ( |3.1| ) and ( ^.3| ), the following relations hold 

L/2 _ ,2t2 



h" 

E' + E^ 
Wh' 
W' + EW 



W 

CiX^h^ + CsXh + C2 
2CiX^h + C3A. 



Substituting the above relations (|3.6| ) for Eq. ( p.5|) , we obtain, 

1)^2, 



■"^riin— 1 



'-'7i,n+l 



1 

12 

^(n 



(3.5) 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 
(3.6e) 



(3.7a) 



(AT- l)(A/'-2) + 6(n- l)(n-Ar) A^ ± -(A/" - 2n + 1)C3A, (3.7b) 



Ar)CiA2 



(3.7c) 



All the other matrix elements are zero. 

The special choices X = ig, Ci = l/2ig, C2 



W{q) 



1 



-sm{gq), E{q) 
9 



l/2ig and C3 



lead to. 



(3.8) 



and correspond to the periodic potential in Ref. 0. We note that Eq. (|3.^ ) is incorporated 
with the perturbation theory defined by Eq. ( p.21| ). We will later carry out non-perturbative 
analysis of this special case in section 



B. Sextic Oscillator Potentials 



Next, we will consider the case where. 



1 



(3.9) 
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This is also a solution of Eq. ( p.9c|) . This special case corresponds to (one of) the cubic type 
in Ref. H]. The Hamiltonians and the supercharge are given by, 



N 



(^f-l)/2 

n 

k=-{j\r-i)/2 



D 



k 



Q-Qo 



(3.10) 
(3.11) 



with, 



W{q)=C,{q-qof + C2{q-qo) + 



C, 



(3.12) 



We note that the Hamiltonians ( 3.10 ) are parity symmetric. The function h{q) can be chosen 

as, 



h{q) 



Bases of the solvable subspaces Vtf are calculated as 



(3.13) 



(3.14) 



Thus, either 0"'" or is normalizable unless Ci = C2 = and the corresponding system Hj^ 
or Hj^ is quasi-exactly solvable. The correspondence between quasi-exact solvability and 
A/'-fold supersymmetry in the case of the sextic potential Eq. (|3.10| ) was recently pointed 



out in Ref. 23 



The non-zero matrix elements of can be obtained by Eq. 
and (|3.13| ), the following relations hold: 



h" 
Wh' 



2h, 

1, 
0, 

4Cih^ + 2C2h + C3, 
W' + EW = 8Cih + 2C2. 

Substituting the above relations ( |3.15| ) for Eq. ( ^.51 ) we obtain, 

S^.-i = 1(^1-1) -2n + 2)T 2C3 

S±„ = ±(Ar-2n + l)C2, 
St.+i = T4(n-Ar)Ci. 

All the other matrix elements are zero. 
If we rewrite Eq. ( |3.12| ) as. 



I). FromEqs. (|3j), {^J^ 

(3.15a) 
(3.15b) 
(3.15c) 
(3.15d) 
(3.15e) 



(3.16a) 

(3.16b) 
(3.16c) 



W{q) = w{q) + 



«^(g) = Ci{q-qo)^ + C2{q-qo), 



(3.17) 



8 



the potential parts Vj^{q) of the Hamiltonians ( p.lO|) are, in terms of w{q), 



(2^3 T A/"- l)(2C3TAr+l) ^ /A/"^ a 



8(g - qo? 



± 



± 



w'{q) + -C2Cs. (3.18) 



We note that in the cases when C3 = {Af ± l)/2 and —{J\f± l)/2, one of the potential-pair 
y^(g) becomes a genuine sixth order polynomial: 



, ,2 4Af±l . 
w{q) + w (g) 

D 

T/-/ A - f ^2 47V± 1 , 

^A^W) = o«^W) ^ w{q) 



C3 
C3 



2 

A/'±l 



(3.19a) 
(3.19b) 



where irrelevant constant terms are omitted. Conversely, a sextic anharmonic oscillator or 
a triple-well potential (with parity symmetry) can be one of the Type A A/'-fold supersym- 
metric pair whenever the potential can be put in one of the form of Eq. ( p.l9| ). When 
C3 = {Af±l)/2, the bases Eq. ( CT ) for V^{q) read. 



—(9 -go) 



2±2exp 



C3 = ^|. (3.20) 



It is worth noting that the solvable subspace consists of the states with definite parity 
(odd for C3 = {Af + l)/2 and even for C3 = {Af — l)/2). We will later see close relation 
between this fact and pattern of the non-perturbative spectral shifts. When C3 = —{Afi: 



l)/2, the bases Eq. ( |3.14|) for Vj^{q) are similarly 



2=^2 exp 



Cl I n4 C*2 / n2 



C3 



A/'± 1 



.(3.21) 



Again, the subspace Vj^ contains only odd-parity states for C3 = — (A/" + l)/2 and only 
even-parity states for C3 = —{Af — l)/2. 

C. Quartic Oscillator Potentials 



In the next, we will consider the case where E{q) = 0. This is also a trivial solution of 
Eq. {^^. From Eq. {^^) we yield. 



W{q) = Ciq^ + C2q + C3. 
In this case, the Hamiltonians and the supercharge are given by 

H^M=lp' + lw{qf±^W'{q), Pu = D 



The function h{q) reads. 



h{q) = q. 



Bases of the solvable subspaces are calculated as. 



(g) = g" exp 



(3.22) 



(3.23) 



(3.24) 



(3.25) 
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Thus, both of 0^ are not normahzable and therefore the system is quasi-perturbatively 
solvable as far as Ci is a non-zero real number. The relation between quasi-perturbative 
solvability and A/'-fold supersymmetry in a special case of the models was pointed out in 
Ref. 0. 

The non-zero matrix elements of can be obtained by Eq. ( ^.51 ). From Eqs. ( |3.22| ) and 
( ^.24| ), the following relations hold: 

h'^ = 1, (3.26a) 
Wh' = Cih^ + C2h + C3, (3.26b) 
W' + EW = 2Cih + C2. (3.26c) 

Substituting the above relations ( |3.26| ) for Eq. ( ^.5[ ) we obtain, 



Si^^_2 = --{n-l){n-2), (3.27a) 

St„-i = Tin - l)Cs, (3.27b) 

Sj„ = ±i(Ar-2n + l)C2, (3.27c) 

S;.+i = T(n-Ar)Ci. (3.27d) 
All the other matrix elements are zero. 



IV. ANALYSIS OF A PERIODIC POTENTIAL 
A. The Valley Method 

Before proceeding to show the results of the analyses, we briefly review the valley 



method Q |3^, |3^, 0, ^ which is employed in this research. For more 

details about the method, see Ref. 0. 

The main problem in quantum theories concerns with the evaluation of the Euclidean 
partition function: 

-S[q] 



Z = Jj Pge-^'^l. (4.1) 

Since the evaluation cannot be done exactly in general, one must find out a proper method 
which enables one to get a good estimation of the quantity. The semi-classical approximation 
is known to be one of the most established methods. Especially, the uses of instantons have 
been succeeded in analyzing non-perturbative aspects of various quantum systems which 
have degenerate vacua |^0|. However, validity of the approximation comes into question 
when the fluctuations around the classical configuration contain a negative mode. Let us 
consider an asymmetric double-well potential as a typical example. For this potential, there 
is a so-called bounce solution as the classical solution which has a negative mode in the 
fluctuations. The negative mode contributes non-zero imaginary part of the spectra in the 
approximation, showing instability of the system. Since the spectra of the model must be 
real, the instability in the approximation must be fake. 

The appearance of a negative mode indicates that the classical action does not give the 
minimum but rather a saddle point in the functional space. In this case, one may expect 
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that the quantity ( |4.1| ) is dominated by the configurations along the negative mode, which 
may intuitively constitute a valley in the functional space. The valley method is a natural 
realization of this consideration. 

At first, we give a geometrical definition of the valley in the functional space g(r) [p^]: 
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5q{r) 



0. (4.2) 



The above definition (|4.2| ) can be interpreted as follows; for each fixed "height" 5'[g], the 
valley is defined at the point where the norm of the gradient vector becomes extremal. 
Introducing an auxiliary field F{q), we can make the valley equation ( [4.2| ) a more perspicuous 
form: 

= Fir), (4.3a) 

jdr' D{t,t')F{t') = AF(r). (4.3b) 
where the operator D is defined as, 

™,,0 = ^!%L (4 4) 

It is now evident that any solution of the equation of motion is also a solution of the valley 
equation (|4^ ) with F(r) = 0. 

Next, we separate the integration along the valley line from the whole functional integra- 
tion. We parametrize the valley line by a parameter R and denote the valley configuration 
by qRir). We then define Faddeev- Popov determinant A by the following: 



dR5 [^j drvR{T)GR{T)J /\\vr] = 1, (4.5) 

where ^r{t) = qij) — qRij) is the fiuctuation over which we will be doing Gaussian inte- 
grations, and Gr{t) is the normalized gradient vector. 



Inserting Eq. ([4.5|) into the functional integral (|4.1|) , expanding the action S[q\ around 
'^nij) = and integrating up to the second order term in ^r{t), we finally obtain the 
one- loop order result: 



Z = J jdR jvqsi^j dT^R{T)GR{T)^ A[^R]e- 



■S[q\ 



where the Jacobian A[y9R] is given by, in this approximation. 
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In the above, det' denotes the determinant in the functional subspace which is perpendicular 
to the gradient vector Gr{t). The valley equation ( |4.3| ) ensure that the subspace does not 
contain the eigenvector of the eigenvalue A. Therefore, we can safely perform the Gaussian 
integrations even when we encounter a non-positive mode. The extension to the multi- 
dimensional valley, which will be needed when there are multiple non-positive eigenvalues, 
is straightforward. 

In this article, we only deal with one- dimensional quantum mechanics where the Euclidean 
action is given by. 



S[q] 



dr 



1 ( dq 



2\dT 



+ V{q) 



(4.9) 



In this case, the valley equations ( [4. 3] ) are explicitly written as, 

d\{T) 



dr^ 



+ V'{q) = F{t) 



d'' 



dr^ 



+ V'\q) 



Fir) 



\F{t). 



(4.10a) 
(4.10b) 



B. Valley-Instantons 

At first, we will analyze a periodic potential. The form of the potential to be analyzed is 
the following: 



1 -,2, ^ e 



V{(1] e) = ^ sin {gq) + - cos{gq). 



(4.11) 



This is a periodic potential with periodicity 2TT/g (unless e = 0) and has two local minima 
at g = 2k7i/g and q = {2k + l)/g (A; = 0, ±1, ±2, . . . ) in one period, see Fig. |I|. 



i 








2 


\ 71 / 





s 





q 



FIG. 1: The form of the periodic potential investigated in this section. 
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Comparing this potential with Eqs. ( |3.2| ) and ( |3.8|) , we find that the system has Type A 
A/'-fold supersymmetry when 

e = ±Ar. (4.12) 

Since the system is defined on a bounded region, all the bases of the solvable space are 
normalizable. Thus, certain linear combinations of them serve as physical eigenstates of 
the Hamiltonian and A/'-fold supersymmetry is not broken dynamically. Therefore, we may 
expect that the non-perturbative corrections for certain A/" physical states will vanish and 
the perturbation series for the corresponding spectra will be convergent when e = iA/". 
We note that the potential ( [4.11|) has the following symmetry: 

v(^q-^;e^ =Viq;-e). (4.13) 

Therefore, we can restrict e to be positive without any loss of generality. 

In the case of e = 0, there are (anti-)instanton solutions of the equation of motion which 
describe the quantum tunneling between the neighboring vacua. The instanton and anti- 
instanton which connect the two vacua at q = kir/g and q = {k + ^)t^ / g are given by. 



qn\T — Tq) = h - arccosf — tanh(r — ro)V (4.14a) 

gQ^^(r — To) = 1 — arccos^tanh(r — tq)) . (4.14b) 



When e 7^ 0, the classical solutions drastically change into the so-called bounce solutions 
which cause fake instability. On the other hand, the solutions of the valley equation (|4.10 



contain a continuously deformed (anti-)instanton which connects the two non-degenerate 
local minima and is called (anti-)valley-instanton 0. 

The solutions of the valley equation ( |4.10|) also contain a family of the configurations, 
which tends to the trivial vacuum configuration in the one limit and tends to well-separated 
valley- instanton and anti- valley- instanton configuration in the other limit. The latter con- 
figuration is called //-valley. The bounce solution is also realized as an intermediate con- 
figuration of this family, which is consistent with the fact that the solution of the equation 
of motion is also a solution of the valley equations. For details, see the numerical result 
in Ref. 0. For the //-valley configuration, it turns out that |A| -C 1 and thus the asymp- 
totic form of the configuration can be obtained by solving the valley equation ( [4.10[ ) with 
perturbative expansion in A: 

g(r) = go(r) + Agi(r) + ■ ■ ■ , F(r) = AFi(r) + X^F^ir) + ■ • ■ . (4.15) 

Indeed, if we denote the distance between the valley-instant on and the anti-valley-instanton 
as /?, the lambda is order A ~ 0(e~^) quantity. The action of the //-valley with the 
boundary condition g(±T/2) = 2kTi / g (T 3> 1) is finally obtained as, 

S^''\R) = S^''\R) = 25j') -^R+ ^(T - /?) - ^e-^ + 0(6"^^), (4.16) 

while the one with g(±T/2) = {2k + l)n/g (T > 1) is, 

S^''\R) = S'-''\R) = 2Sl!^ + -R--(T-R)- ^e-"^ + 0(e-'^), (4.17) 

2 2 g^ 
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where Sq denotes the Euchdean action of one (anti-)instanton Eq. ( 4.14| ) and amounts to, 



= 4- (4-18) 
9 

In Eqs. ( [4.16| ) and ( |4.17| ), the fourth term can be interpreted as the interaction term between 
the valley-instanton and the anti-valley-instanton. Therefore, the minus sign indicates that 
the interaction is attractive. 

The other type of the solutions emerges in this case, which is asymptotically composed 
of two successive valley-instantons or two successive anti-valley-instantons. We call them 
//-valley and //-valley, respectively. These configurations do not appear in the case of 
double-well potentials since they connect every other vacuum. The Euclidean action of 
them with large separation R can be also calculated in the same way as, 

S(''\R) = S^''\R) = 2S'S') -^R+ Ut -R) + -^e~^ + 0{e-^^), (4.19) 

Z Z g 

for the configuration with g(-T/2) = 2kTT/g and g(T/2) = {2k ± 2)7v/g (T > 1), and, 

S^''\R) = S^"\R) = 2^J') + ^R- Ut -R) + 4e-^ + ©(e-^^), (4.20) 

Z Z g 

for the one with q{-T/2) = {2k + l)Tr/g and q{T/2) = {2k + l±2)'K/g (T > 1). Note that the 
sign of the fourth term is plus and thus the interaction between the (anti-)valley-instantons 
in this case is repulsive. 



C. Analysis of Two-valley Sector 

From the results on the //-valley, the contribution of the //-valley to the partition func- 
tion can be written as the following form: 



4g-T/2 rT 

IR ('i ' - J 

dtj^{t)e~'^'\ (4.21) 



Jo 

4e-r/2 



where we have changed the integration variable R to t = g^S^\R) in the second line. The 
integration contour Cy is [0, 2g'^SQ^) and the integrand jF(t) has a singularity at t = 2(7^5'q^\ 
The integral Eq. ( [4.21| ) contains both the perturbative contribution at t ~ and the non- 



perturbative one at t ~ 2g'^SQ\ To separate the perturbative and the non-perturbative 
contribution, we deform the contour Cy to the sum of Cp and Cnp: 



-^9 Jcp JCnp 

= Zf\g^) + Zlll\g% (4.22) 



as is shown in Fig. ^. 
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FIG. 2: Deformation of the contour Cy to the sum of Cp and Cnp- 



We identify the first term as the formal Borel summation of the perturbation series and the 
second term as the non-perturbative contribution. For the non-perturbative contribution, 
the following analytic property holds. If we perform the analytic continuation of Z-^p{\g'^\e^^) 
from 9 = Oto 9 = n, the contour for Znp changes from Cnp(O) to CNp(vr), as shown in Fig. |^. 
In the weak coupling limit, the integral of CNp(7r) can be well-approximated by that of Cy 
because the dominant contribution of the integral comes from t ~ 2g'^SQ\ Therefore, in the 
case of g'^ = \g'^\e^'" when the interaction between valley- instantons is repulsive, the following 
relation holds approximately: 

Z^p{\g'\en ^ Z{\g'\e'^). (4.23) 

This relation coincides with what Bogomolny suggested heuristically as a method of evalu- 
ation of the instanton-anti-instanton contribution [HT| . 




FIG. 3: The change of the contour Cnp(^) as 9 is changed from zero to vr. 



An immediate consequence of our decomposition of the perturbative and non-perturbative 



15 



contribution is, 



ImZp + Im^NP = 0, (4.24) 
since Z = Zp + Znp is real. From the relation above, the dispersion relation becomes 0, 

1 / , Zp(z) 



if/'f-^^. (4.25) 



where Cg2 is the counter around z = g"^ and we have neglected the contribution from the 
singularities far from the origin. The last line of Eq. (f4.25|) gives the relation between 



the coefficients of the perturbation series and the imaginary part of the non-perturbative 
contribution. Rewriting the above relation in terms of the energy spectra, we find for the 
perturbative part of the spectra Ep{g'^) = YlT=o^^^^3'^^ i ^^^^ coefficients a^'') can be 
estimated as. 



The situation in the case of //-valley is completely the same as that in the //-valley. 

On the other hand, the situation in the case of //-valley is different, reflecting the fact 
that the //-valley conflguration cannot be deformed continuously to the trivial vacuum 
conflguration. The contribution of the //-valley to the partition function has quite similar 
form to that of the //-valley: 

Z^'') = ^^— dt:F{t)e~"3\ (4.27) 

However, the integration contour is now Cv = {2g'^S'i\g'^S^^^\'d)] and thus is disconnected 
to the perturbative region t ~ 0, see Fig. ^ 

This means that Eq. ( [4.27 ) contains only the non-perturbative contribution. 



Z^g') = Z{g'). (4.28) 

Therefore, we need not separate the integration as in the case of the //-valley. As a con- 
sequence, the //-valley conflguration does not contribute to the imaginary part and also to 
the large order behavior of the perturbation series, which will be conflrmed in the examples 
in sections |^ and 0. Furthermore, since the interaction for the //-valley conflguration is 



repulsive, as has been observed in Eqs. ( 4.19|) and ([4.20 ), the integral is dominated around 



can be approximated by the integral on the contour \2g'^SQ \ oo) in the weak 
coupling limit. Therefore, analytic continuation is not needed. The situation in the case of 
the //-valley is completely the same as that in the //-valley. 



D. Multi-valley Calculus 



Utilizing the knowledge of the (anti-)valley-instantons and the interactions between them 
obtained previously and applying the manipulation discussed in subsection [IV C|, we will 
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2oW 



2gX 



8Y'(0) 



FIG. 4: The integration contour Cy in the case of the //-valley. 



evaluate the partition function Z = tie'^^ by summing over those configurations made 
of several (anti-)valley-instantons which satisfy a boundary condition in T. The periodic 
boundary condition for a configuration g(r) is in general given by 

q{T + T)=q{r). (4.29) 



For a system which has a periodic potential like Eq. ( |4.11| ), however, the condition ( |4.29| ) 
can be relaxed and be replaced with, 

2k7T 

q{T + T)=q{T) + (A; = 0,±l,±2,...). (4.30) 

9 

This condition restricts the number of the valley-instantons to be even, 2n. The non- 
perturbative contributions from the multi-valley configurations satisfying Eq. ( [4.30| ) can be 
calculated by the extension of the technique developed in Ref. [|^. We divide the time 
interval < r < T into n regions and put a valley-instanton pair on each of the region. In 
order to distinguish what kind of pairs, we introduce the indices and q for the z-th region 
(mod n) as follows: 

(i) ej = 1, ii = 1 for //-valley, 

(ii) ej = 1, ii = —1 for //-valley, 

(iii) ej = —1, ii = 1 for //-valley, 

(iv) ej = — 1, ii = —1 for //-valley. 

In this way, the allowed configurations for a given valley-instanton number 2n are exhausted 
by the allowed combinations of the set {e^, Cj} (i = 1, . . . , n). Combining the results on well- 
separated valley-instanton pairs with the above conventions, the well-separated multiple 
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valley-instanton action for given n and {ej, ej} is expressed as, 



9 



i=l 



9' 



i=l 



(4.31) 



i=l 



i=l 



where Ri is the distance between the {2i — l)-th and 22-th (anti-)valley-instanton and Ri the 
one between the 2z-th and the {2i + l)-th (anti-) valley-instanton (mod n) and e^+i = ei, see 
Fig. |. 



2m+l 




2m- 1 



2m-2 



FIG. 5: The collective coordinates Ri and Ri for a 2n valley-instantons configuration. 



The sum of the contributions from the 2n valley-instantons configuration can be written 

as, 



vNP — 

n=l 



(4.32) 



where denotes the contribution of the Jacobian and the i?-independent part of the de- 
terminant for one valley-instanton-pair and is calculated as, in this case, 



a' = —,e-^l'\ (4.33) 
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The term J„ is given by, 

r>00 / " 



J = ^ 

n 







,i=l 



vi=l 



1 + e" 1— e" 



1=1 



i=l 



eieiC e^ei+ie - i- ^(e^ + e^) ) . (4.34) 

i=l ^ i=l j=l 

To calculate the sum over the set {ej,ej}, we introduce the following transfer matrices: 



exp( 



_8_p — fii 



I'd) exp(4e 



8 ^-Rr 



g ' ^9 



8 



exp(-^e ^0 exp(^e" 
exp(^e-^0 exp(-^e--^0 



(4.35a) 
(4.35b) 



Then, using these matrices we have 

pico — Tj 

— / dse"""tr I I I 

-Ts 



Jn 



271171 

T 



27rm 



■loo—rj 
ioo—rj 



rfse-'^tr [T(s,e,79)' 



tOO — f] 



where, 

K(s - e/2)e-'^'' /(s - e/2) 
/(s - e/2) K{s - e/2)e'^ 

In the above, K and / are defined by. 



r(s,e,^9) 



K(s + e/2) /(s + e/2) 
/(s + e/2) ir(s + e/2) 



where the manipulation explained in subsection [iV C| is utilized for estimating each of the 
integration. The calculation of the trace in Eq. ( |4.36| ) can be done by diagonalizing T. If 
we denote the two eigenvalues of T as t±, we immediately yield. 



(4.36) 
(4.37) 

(4.38a) 
(4.38b) 



J n 



T 



27rm 



/ZOO — ?y 
c/se-^^(t+(s)" + t_(s)"V 
'ioO—Tj 



(4.39) 



From Eq. ( [4.37|) , t± is evaluated as, 

t±(s) = K{s + e/2)K{s - e/2) cosd + I{s + e/2)I{s - e/2) 



± 



+ 



I{s + e/2)' - K{s + e/2)' K{s - e/2)' sin' ^9 



K{s + e/2)I{s - e/2) + I{s + e/2)K{s - e/2) cosd 



1/2 



(4.40) 
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Finally, combining Eqs. ( [4 .321) and (|4.39| ) we obtain the non-perturbative contribution to 
the partition function: 



T 

27ri 



dse-^' ln( 1 - a'HJs) ) ( 1 - a'HJs 



—loo—rj 



(4.41) 



E. Non-perturbative Contributions 



From the results in Eqs. ( [4.38| ), ( [4.40|) and (|4.41| ), the non-perturbative contributions to 
the spectra are determined by the following equation: 



SA"*-*" / 1 6 



2 2 



(4.42) 



where, 



f3±{E,e,^) = cos'd + 



(i^-l-f )2 + (_)(i?-f + f ) 



m + (_)(i?-^)22 cos^ - sin^ (4.43) 
We will solve the above equation by the series expansion in a: 

+ E^}^ + ■ ■ ■ 



n± 



(4.44) 



where i?„_^ stands for the spectra corresponding to, in the limit (7 ^ 0, the eigenfunctions 
of the shallower potential wells and En_ for the ones corresponding to the eigenfunctions of 
the deeper potential wells. For e ^ M {M = 1, 2, 3, ... ), the 1st order contributions vanish 
and the leading 2nd order contributions are calculated as follows: 



Eg = 2[cost9 + (-)±^] 



\n±+l 



2n±±€ 



r(-n±Te) 



(4.45) 



For e = M {M = 1, 2, 3, ... ), all the harmonic spectra En} of the shallower wells and the 

higher harmonic spectra Enl of the deeper wells degenerate for n_ = n+ + A/", see Fig. |^ 
for e = A/" = 1, 2. Between these degenerate states, resonant tunneling enhances the non- 
perturbative corrections, and results in order contributions to the spectra En_ and i?„_^ 



with n_ 



+Ar: 



n± 



^(2) 

n± 



2[l + (-)^cos?9] f 8 



(1)2 



n± 



21n(- 



1 + (-)^ COS^ 



g- 

+ 2 In 



9' 



(4.46) 
(4.47) 



where il^{z) = d\nr{z)/dz is the digamma function. For the other spectra, say, the lower 
En_ with ?T,_ < A/", the contributions are the same as Eq. ( [4.45|) . We see from Eq. ( |4.45|) 
the non-perturbative corrections for these lower En_ vanish, at least up to order a^, at 
■(? = (tt) when e = is odd(even) and thus, from Eq. ([4.12| ), the system has A/"- fold 
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super symmetry. This means that when the system is TV-fold supersymmetric with odd(even) 
TV, among the physical states for each of the lower TV spectral bands, the state which satisfies 
the periodic (anti-periodic) boundary condition does not receive non-perturbative correction. 
From Eq. ( p^ ) in this case, these physical states are surely the elements of the solvable 
subspace V^. Therefore, the results are consistent with the fact that TV- fold supersymmetry 
in this case is not broken dynamically. 




FIG. 6: Degeneracies of the harmonic spectra for (a)e = 1 and (b)e = 2. 



Finally, we make a remark on the resultant equation ( ^4.42|) . The origin of the disap- 
pearance of the non-perturbative corrections discussed above comes from the factor (3± in 
Eq. (|44^) . Indeed, to make the l.h.s. of Eq. (|442|) finite when (3± = 0, the gamma functions 
must diverge adequately. This happens only when E = En± for certain values of n±, and 
therefore the non-perturbative corrections must vanish.^ From the derivation of Eq. ( [4. 421) , 
we can see that the appearance of the factor f3± is achieved by taking into account both 
of the repulsive and attractive interactions between the valley-instantons properly. There- 
fore, we guess that naive application of the dilute-gas approximation can hardly lead to the 
correct results even for the ground-state energy. 



F. Large Order Behavior of the Perturbation Series 

The large order behavior of the perturbation series in for the spectra can be estimated 
by the same way as in the case of the double-well potential. From the non-perturbative 
contributions Eqs. ( |4.45| )-( |4.47| ), we can easily see that the imaginary parts of them are 
continuous, at least up to the order a^, in e and yield. 



27r 



2n±±e 



n±\ T{n± + 1 ± e) \g 

which are valid for arbitrary e. Then, if we expand the spectra in power of g"^ such that 



(4.48) 



(4.49) 



r=l 



Note that the argument does not depend on the perturbative expansion in a. 
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the large order behavior of the coefficients a^'^^ for sufficiently large r are calculated as, using 
Eqs. and iWM) . 



a2~^n±(e) 4-'T(r + 2n± + l±e) 



where, 



^"^^"^ 7rn±!r(n± + l±e)- 



(4.50) 



(4.51) 



Equation ( |4.5(J| ) shows that the perturbative coefficients diverge factorially unless the pref- 
actor A{t) vanish. From Eq. ([4.51|) , we can find the disappearance of the leading divergence 
takes place only when e = Af {Af = 1, 2, 3, ... ). Comparing the results with Eq. (|4.12|) and 
taking the symmetry (^4.13|) into account, we see that the above cases completely coincide 
with the case where the system possesses Type A AA-fold supersymmetry. Therefore, the 
results of the valley method analyses are consistent with a consequence of Type A AA-fold 
supersymmetry, that is, the non-renormalization theorem. 

V. ANALYSIS OF A TRIPLE- WELL POTENTIAL 

In this section, we will analyze a sextic triple-well potential. The form of the potential 
to be analyzed is the following: 



This has three local minima at g = and g ^ ±1/(7 for eg"^ <^ 1, see Fig. 0. 



(5.1) 




FIG. 7: The form of the triple-well potential investigated in this section. 



Comparing this potential with Eq. ( |3.19| ), Ci, C2 and go in w{q) being, 

Ci = -9^ C2 = 1, go = 0, 
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(5.2) 



we can easily see that the system has Type A TV-fold supersymmetry when 

AU±l , ^ 

e = ±^— . (5.3) 

More precisely, the system becomes one of the A/'-fold supersymmetric pair H^; when 
e = {AN ± l)/3 and Hjj- when e = —{AN ± l)/3. As has been explained in subsection 
111 B| , A/'-fold supersymmetry does not break in the cubic type because either bases of the 



solvable subspace Vjf or those of Vj^ are normalizable in general. Since Ci < in this 
case, bases of the solvable subspace are normalizable and physical while those of are 
not. Therefore, we may expect that the non-perturbative corrections for certain N states 
will vanish when e = {AN ± l)/3 while those will not when e = —{AN ± l)/3, though the 
perturbation series for the corresponding spectra will be in both the cases convergent. 
As far as we know, little has been investigated for triple-well potentials. We have found 



only two references [^, Q on the subject, both of which employed the dilute-gas approxima- 
tion. However, from the consideration mentioned at the end of subsection |1V E| , we intend to 
analyze beyond the dilute-gas approximation using the same technique of the valley method 
as in section P^ .^ 

In the case of e = 0, the three local minima of the potential have the same potential 
value. Thus, there are (anti-)instanton solutions of the equation of motion which describe 
the quantum tunneling between the neighboring vacua: 

(/■) 1 1 \ 1 

^0 ^^~^°^ = ^^ (l + eT2(^-o))i/2 ' ^0 (^-^o) = ^- (i + e^^(— o))i/2 - (^-^^ 

When e 7^ 0, the solutions of the valley equation now become the (anti-)valley-instantons. In 
this case, there are three kinds of the solutions of the valley equation which are asymptoti- 
cally composed of two (anti-)valley-instantons. Contrary to the periodic potential in section 
[rV| , there are two different //-valley or //-valley configurations in this case since the curva- 
ture at the central potential bottom (at g = 0) is different, even at the leading order of g"^, 
from the one at the side potential bottoms (at q ~ ±1/(7); the // (//)-valley which satisfy 
g(±T/2) = (T ^ 1) are different from the ones which satisfy g(±T/2) ~ 1/g or —1/g 
(T 3> 1). The Euclidean action of the former with large separation R can be calculated by 
the perturbative expansion in A ~ 0{e~^^) as follows: 

S^''\R) = S^''\R) = 2Sl'^ -eR+UT - R)- ^e"'"^ + 0(e-^^), (5.5) 

2 g 



while the one of the latter with large separation R can be calculated in the same way as, 
S^"\R) = S^''\R) = 25j'^ + ^R- e{T - R) - ^e'"^ + 0(6^'^), (5.6) 

where 5*0^^ denotes the Euclidean action of one (anti-)instanton Eq. ( |5.4D and amounts to. 



Si>" = ^- (5.7) 



^ Indeed, the results of Ref. Q for the lowest three energies are different from ours and are not consistent 
with A/'-fold supersymmetry. 
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The other type is //-valley or //-valley. The Euclidean action of them with large separation 
R can be also calculated in the same way as, 

S^''\R) = S^''\R) = 2Si'^ + -R- e(T - R) + ^e""^ + 0{e~^^). (5.8) 

2 

A. Multi-valley Calculus 

The evaluation of the partition function Z = tre~^"^ by summing over multi-valley- 
instanton configurations can be done in the same manner as those for the double-well and 
the periodic potentials. One can easily see that in order to incorporate with the periodic 
boundary condition in T, the number of the valley-instantons in a period T must be even. For 
a given number 2n of the valley-instantons, however, there are still several configurations. 
If we regard a configuration as n valley-instanton pairs, we have four kinds of pair, //-, 
//-, //-, and //-valleys. We denote the number of the //- and //-valley as riu and ri//, 
respectively, and that of the others as rijj. Contrary to the periodic potential case, the 
particle must come back to the start point after the period T passes in this case. Therefore, 
we must impose Eq. ( |4.29| ) rather than Eq. ( [4.301 ). This condition results in ri// = rijj. As 
a consequence we have, 

2n// + Ujj = n. (5.9) 

This restriction shows that for a given n there are [n/2] + 1 variety of the riu value. For n 
and nil fixed, however, the configuration is not determined uniquely yet. There remains a 
freedom of the permutation of the pairs. The number of cases can be calculated if one notices 
that the configuration is uniquely determined as far as the position of the //- and //-valleys 
among n area is fixed. We denote a set of the position as {in}- It is therefore clear that for 
given n and njj there are nC2nji configurations of the multiple valley-instantons. Combining 
the results on well-separated valley-instanton pairs with the above considerations, the well- 
separated multi- valley-instanton action for given n, nn and {in} is expressed as, 

n n 

e 



2 

i=l 1=1 

^ i=i ^ ^ mill} 

where Ri is the distance between the {2i — l)-th and 22-th (anti-)valley-instanton and Ri the 
one between the 2i-th and the {2i + l)-th (anti-) valley-instanton (mod n), see Fig. |^. 

The sum of the contributions from the 2n valley-instantons configuration can be written 

as, 

oo [n/2] 

Znp = ^ a^V^, Jn=Y.Y. (5-11) 

n=l nii=0{iij} 

where denotes the contribution of the Jacobian and the /^-independent part of the de- 
terminant for one valley-instanton-pair and is calculated as, in this case, 

= yle-'/'^\ (5.12) 
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The term Jnln\] is given by, 



q{iii} — 
n,nii 



T 



n 



POD / \ / _\ 



exp 



, i=l 



i=l 



i=l 



ie{j//} 



(5.13) 



In the above expression, we notice that for n and nn fixed, the contribution J^^^nn does not 

(5.14) 



depend on the choice of the set {in}. This means the following equality, 

{in} 



/ J '-'n,nii 
{in} 



n 
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where Jn.nu is the contribution Jn%u for a specific {in} and is evaluated as, 

rj-f poo 



J. 



n,njj 



{in} 

i=i J \i=i 
6 [ + exp 



n 



j=i 



i=l ^ ^ i=\ 



-Ri 



3 .=1 

rp I'ioo—rj 



27rm 



dse-^'K_{sYK^^\sf''"Kf{s 



1=1 i=2nii+l 



(5.15) 



In the last expression for Jn,nn^ several i^''s are defined by, 



Kjs) 



dR exp 







(s - 1 + e)i? + — e 
9^ 



'2R 



1 / 1 ^"5+5 / s 1 e 



(5.16a) 



K^^\s) = I dRexp 



1 ~ 2 . 

2 2 



2^22 



dR exp 



r(-. + i + i 

' 2 2 



s - - - - R+-,e-^ 



(5.16b) 



2 2 



2 2/ 1 e 
' 2 2 



where the manipulation explained in subsection p.V C| is again utilized for estimating 
the integrations. Eventually, from Eqs. ( [5.11|) , (|5.14| ) and ( p.l5| ) we obtain, 

[n/2] 



(5.16c) 
each of 



f n \ 

Znp = I] (2nfr)^™ 

n=l nii=0 ^ ^ 

'dse-^Mn (l - a^K4s)K^^\s)^ (l - a2ir_(s)iri"^(s)) 



loo—rj 



where. 



(5.17) 
(5.18) 



B. Non-perturbative Contributions 

From the results in Eqs. (|5.16|) - (|5.18|) , the non-perturbative contributions to the spectra 
are determined by the following equation: 

2\^"2-2 / 1 e\ / 1\ 2-5+2 r E 1 e 



a^P.(E,e){ j,j r(^-£ + ^ + -j(^~l r I -| + ^ - - j = 1, (5.19) 
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where, 



P±{E,e)=^-^^ . (5.20) 

We will solve the above equation by the series expansion in a: 

= + aE« + a'Ei^J + ..., = no + ^ + |, (5.21a) 

E^^ = + aEl^l + a'Ei^l + ■■■, = 2n± + 1 - e, (5.21b) 

where Eng stands for the spectra corresponding to, in the limit g ^ 0, the eigenfunctions 
of the center potential well and En^ for the ones corresponding to the parity eigenstates 
obtained by the linear combinations of the eigenfunctions of the each side potential well. 
For e 7^ ±{2J\f + l)/3 (A/" = 0,1,2,...), the 1st order contributions vanish and the leading 
2nd order contributions are calculated as follows; 




(5.22) 
(5.23) 



In this case, degeneracies of the harmonic oscillator spectra for the each potential well only 
occur between the both side wells. The different non-perturbative contributions for n± in 
Eq. ( |5.23| ) show the splitting of the degeneracies via the quantum tunneling as in the case 
of symmetric double-well potentials. 

When e = {4JV + l)/3 (A/" = 0,1,2,...), all the even-parity central harmonic spectra 
E^^ and the higher side harmonic spectra En} degenerate for n± = itlq + A/", see Fig. |^(a) 
for e = 5/3 {U =1). 




FIG. 9: Degeneracies of the harmonic spectra for (a)e = 5/3 and (b)e = —1. 

It is interesting, however, that the interference due to the quantum tunneling only occurs 
between the same (even-)parity states. As a consequence, E2mo and Enj^ satisfying n+ = 
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mo + M acquire order contributions as follows: 



Ti+/2mo 



±1 



n+!(2mo)! \g 



2mo 



9' 



(5.24) 



E, 



(2) 

n_|_/2mo 



E 



(1)2 
n-|-/2mo 



4 

+ ln 



In 



2^/'(2mo + 1) 



(5.25) 



and the lower E„ 



For the other spectra, say, £"„_, £"27710+1 "^^^ ^^^^ ±v^vvv.x 
are the same as Eqs. (|]2|) and (|Jl). When e = -(iA/"- 1)/3 (TV 



with n+ < TV, the contributions 
1, 2, 3, ... ), all the side 

harmonic spectra Enl and the higher even-parity central harmonic spectra -E2mo degenerate 
for mo = n± -|-TV, see Fig. ^(b) for e = — 1 (TV =1). In this case, the interference also occurs 
only between the same (even-)parity states. The contributions for E2mo and En+ satisfying 
mo = n+ + TV are given by the same as Eqs. ( |5.24|) and (|5.25|) . For the other spectra, say, 
En^, E2mo+i and the lower E2mo with mo < TV, the contributions are the same as Eqs. (|5.22|) 



and (lOal ). 

When e 

n(0) 



(4TV- l)/3 (TV = 1,2,3, 



all the odd-parity central harmonic spectra 

(0) 



i?2mo+i and the higher side harmonic spectra Enl degenerate for n± = mo+TV, see Fig. p!0|(a) 
for e = 1 (TV= 1). 




FIG. 10: Degeneracies of the harmonic spectra for (a)e = 1 and (b)e = —5/3. 



In this case, only the odd-parity states interfere and yield order contributions for 
n_ = mo + TV: 



E 



(1) 

n_/2mo+l 



±1 



n_.!(2mo + 1)! \g 



2mo+l / \ 



9' 



(5.26) 



E 



(2) 

n_/2mo+l 



^(1)2 

n_ /2mo+l 



+ ln 



In ^ + In 



9"/ \9- 
^{n_ + 1) - 2^(2mo + 2) 



(5.27) 
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The contributions for the other spectra, say, En^ , -E'2mo and the lower En_ with n_ < TV, are 
given by the same as Eqs. (|5:22|) and (lOl) . When e = -{4Af + l)/3 (TV = 0, 1, 2, . . . ), all 

the side harmonic spectra En} and the higher odd-parity central harmonic spectra -E'2m(,+i 
degenerate for mo = n-t + TV, see Fig. |10|(b) for e = —5/3 (TV = 1). Only the odd-parity 
states interfere in the same way and the non-perturbative contributions for -E2mo+i ^ind En_ 
satisfying mo = n_ + Af are the same as Eqs. ( |5.26| ) and ( |5.27| ). Again, the expressions 
for the other spectra, say, -E„_|_, -E2mo and the lower -E2mo+i with mo < TV, are given by 
Eqs. QO^ ) and ( ^ . 

From the whole result obtained here, we see that the non-perturbative corrections vanish 
only when e = {4J\f ± l)/3 (TV = 1,2,3,...). More precisely, when e = (4TV + l)/3, 
Eq. ( |5.23| ) is applied for the even-parity states labeled by the quantum number n_ and 
results in Enl = for all n_ < TV. Similarly, when e = (4TV — l)/3, Eq. ( |5.23| ) is applied 

(2) 

for the odd-parity states labeled by the quantum number and results in En+ = for all 
n+ < TV. It should be noted that in the case of e = -(4TV ± l)/3 (TV = 1, 2, 3, ... ) the 
non-perturbative corrections do remain although the models are TV-fold supersymmetric, 
reflecting the fact that they are only quasi-perturbatively solvable but are not quasi-exactly 
solvable. These results are just what we have expected from the general properties of TV-fold 
sup er sy mmet ry. 



C. Large Order Behavior of the Perturbation Series 

The large order behavior of the perturbation series in for the spectra can be estimated 
by the same way as in the case of the double-well and periodic potentials. From the non- 
perturbative contributions Eqs. ( [5.22| )-( ^.27| ), we can easily see that the imaginary parts of 
them are continuous, at least up to the order a^, in e and yield, 

TT /2\"Vl \ + 

ImEno 5 ^ — — , (5.28a) 

TT / 2 \ ^"±^^"5*^ / 1 \ 

ImE„, — ; — -. 5 ^ ( — ( ^ , (5.28b) 

n±!r(2n± + |-|e) V^7V V^V 

which are valid for arbitrary e. Then, if we expand the spectra in power of such that, 

oo 



r=l 



the large order behavior of the coefficients a'-''^ for sufficiently large r are calculated as, using 
Eqs. dg) and iWM) . 



a^:l ~ An, (e) rT[r+^-no + \ + \^ , (5.30a) 
~ (e) 2^ r f r + 3n± + ^ - , (5.30b) 
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where, 



TT n±! r (2n± + f - f ej 

Equations ( [5 .301 ) show that the perturbative coefficients diverge factorially unless the prefac- 
tor v4(e)'s vanish. From Eq. ( [5.31[ ), we can find the disappearance of the leading divergence 
takes place only when e = ±(2n + l)/3 (n = 1,2,3,...). More precisely, we obtain the 
following results: 

(1) e=(4Ar±l)/3(Ar=l,2,3,...) 

An± (e) = for n± < M . 

(2) e = -(4Ar + l)/3 (AT = 1, 2, 3, ... ) 

^2mo+i(e) = for mo < M. 

(3) e = -(4Ar-l)/3(Ar=l,2,3,...) 

^2mo(e)=0 for mo<A/'. 



Comparing these results with Eq. ( |5.3| ), we see that the above cases completely coincide 
with the case where the system possesses Type A A/'-fold supersymmetry. Again, the results 
of the valley method analyses are consistent with the non-renormalization theorem. 



VI. SUMMARY 



In this article, we have made non-perturbative analyses on the models which can be A/'- 
fold supersymmetric at specific values of the parameter. Combining the results obtained in 
this article with the ones in Ref. p[, we get the following: 

(1) For all the potentials investigated (double- well, triple- well, periodic), the leading di- 
vergence of the perturbation series disappears when and only when they are A/'-fold 
supersymmetric. The non-renormalization theorem ensures that A/'-fold supersymme- 
try is sufficient for the disappearance of the divergence. The results indicate that it 
may also be necessary. 

(2) The non-perturbative corrections to the spectra for certain states vanish when and 
only when the models are quasi-exactly solvable (triple- well, periodic). 

(3) For the quasi-perturbatively solvable potentials (double- well, triple- well), the non- 
perturbative corrections remain although they are A/'-fold supersymmetric. 
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As was mentioned in Ref. the quasi-solvable models constructed by s/(2) generators do 
not always have normalizable solvable states. Although the conditions on the normalizability 
of the models were fully investigated in Ref. it remains unclear what is the role of the 
partial algebraization of the models without normalizable solvable states. The results listed 
above provide an answer to this problem. Even though the solvable wave functions are not 
normalizable, they can be normalizable and thus make sense in the perturbation theory. In 
this case, the spectra corresponding to the solvable states also make sense in the perturbation 
theory. As was shown in Ref. the perturbation series for them are convergent since they 
are the solutions of a finite order algebraic equation. However, the fact that the solvable 
states and the corresponding spectra make sense only in the perturbation theory inevitably 
means the existence of the non-perturbative effects, which is in contrast to the case of the 
quasi-exactly solvable models. That is why we have called the case quasi-perturbatively 
solvable. 

Finally, we would like to mention about applicability of the dilute-gas approximation. As 
has been mentioned previously, the dilute-gas approximation cannot give proper results, that 
is, consistent results with AT-fold supersymmetry, for both the potentials Eqs. ( |4.11|) and 
Therefore, it seems that the success of the dilute-gas approximation for the symmetric 
double-well potential is rather exceptional and apphcability of it is quite limited. 
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